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We theoretically explore coherent information transfer between ultra-slow light pulses and Bose- 
Einstein condensates (BECs) and find that storing light pulses in BECs allows the coherent conden- 
sate dynamics to process optical information. We consider BECs of alkali atoms with a A energy 
level configuration. In this configuration, one laser (the coupling field) can cause a pulse of a second 
pulsed laser (the probe field) to propagate with little attenuation (electromagnetically induced trans- 
parency) at a very slow group velocity (~ 10 m/s), and be spatially compressed to lengths smaller 
than the BEC. These pulses can be fully stopped and later revived by switching the coupling field 
off and on. Here we develop a formalism, applicable in both the weak and strong probe regimes, 
to analyze such experiments and establish several new results: (1) We show that the switching can 
be performed on time scales much faster than the adiabatic time scale for EIT, even in the strong 
probe regime. We also study the behavior of the system changes when this time scale is faster than 
the excited state lifetime. (2) Stopped light pulses write their phase and amplitude information onto 
spatially dependent atomic wavefunctions, resulting in coherent two-component BEC dynamics dur- 
ing long storage times. We investigate examples relevant to Rb-87 experimental parameters and see 
a variety of novel dynamics occur, including interference fringes, gentle breathing excitations, and 
two-component solitons, depending on the relative scattering lengths of the atomic states used and 
the probe to coupling intensity ratio. We find the dynamics when the levels \F — 1, Mf = — 1) and 
\F — 2, Mf — +1) are used could be well suited to designing controlled processing of the information 
input on the probe. (3) Switching the coupling field on after the dynamics writes the evolved BEC 
wavefunctions density and phase features onto a revived probe pulse, which then propagates out. 
We establish equations linking the BEC wavefunction to the resulting output probe pulses in both 
the strong and weak probe regimes. We then identify sources of deviations from these equations 
due to absorption and distortion of the pulses. These deviations result in imperfect fidelity of the 
information transfer from the atoms to the light fields and we calculate this fidelity for Gaussian 
shaped features in the BEC wavefunctions. In the weak probe case, we find the fidelity is affected 
both by absorption of very small length scale features and absorption of features occupying regions 
near the condensate edge. We discuss how to optimize the fidelity using these considerations. In 
the strong probe case, we find that when the oscillator strengths for the two transitions are equal 
the fidelity is not strongly sensitive to the probe strength, while when they are unequal the fidelity 
is worse for stronger probes. Applications to distant communication between BECs, squeezed light 
generation and quantum information are anticipated. 

PACS numbers: 03.67.-a, 03.75.Fi, 42.50 



I. INTRODUCTION 

In discussions of quantum information technology , 
it has been pointed out that the ability to coherently 
transfer information between "flying" and stationary 
qubits will be essential. Atomic samples are good can- 
didates for quantum storage and processing due to their 
long coherence times and large, controllable interactions, 
while photons are the fastest and most robust way to 
transfer information. This implies that methods to trans- 
fer information between atoms and photons will be im- 
portant to the development of this technology. 

Recently the observation of ultra-slow light (USL) 
propagating at group velocities more than seven or- 
ders of magnitude below its vacuum speed (V g ~ 10~ 7 c), 
and the subsequent stopping and storing of light pulses in 



atomic media 

ELI 

has demonstrated a tool to possibly 
accomplish this [6j. The technique relies on the concept 
of electromagnetically induced transparency (EIT) [3 in 
three level A-configuration atoms (see Fig.^a)). A cou- 
pling light field fl c is used to control the propagation 
of a pulse of probe light tt p . The probe propagates at 
a slow group velocity and, as it is doing so, coherently 
imprints it's amplitude and phase on the coherence be- 
tween two stable internal states of the atoms, labelled 
|1) and 1 2) (which are generally particular hyperfine and 
Zeeman sublevels). Switching the coupling field off stops 
the probe pulse and ramps it's intensity to zero, freezing 
the probe's coherent (that is, intensity and phase) infor- 
mation into the atomic media, where it can be stored for 
a controllable time. Switching the coupling field back on 
at a later time writes the information back onto a revived 
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FIG. 1: Schematic of ultra-slow light (USL) and 
stopped experiments, (a) The A energy level structure 
which we consider has two stable states, |1) and |2), and an ex- 
cited state 1 3) which decays at a rate T (which is 10 MHz 
for Na and (2tt) 6 MHz for Rb-87). Atoms which sponta- 
neously decay from |3) are assumed to exit the levels under 
consideration, (b) We consider the two light fields to be co- 
propagating and input in the +z direction (the long axis of the 
Bose-Einstein condensate (BEC)) at Zj n . The output intensity 
at 2 ou t can be detected experimentally with a photo-multiplier 
tube (PMT). The probe field is pulsed with half- width r . The 
coupling beam is c.w. but can be switched off quickly to stop 
the probe pulse while it is in the BEC. Switching it back on 
later regenerates the probe light pulse. 



probe pulse, which then propagates out of the atom cloud 
and can be detected, for example, with a photo-multiplier 
tube (PMT) (see Fig. ^b)). In the original experiment 
0, the revived output pulses were indistinguishable in 
width and amplitude from non-stored USL pulses, indi- 
cating the switching process preserved the information in 
the atomic medium with a high fidelity. 

In addition to the ability to store coherent information, 
compressing ultra-slow probe pulses to lengths shorter 
than the atomic cloud puts the atoms in spatially depen- 
dent superpositions of |1) and |2), offering novel possibil- 
ities to study two component BEC dynamics |8|. In |9j 
slow light pulses propagating in Bose-Einstein conden- 
sates (BECs) [l(| of sodium were used to create density 
defects with length scales ~ 2 /im, near the BEC healing 
length, leading to quantum shock waves and the nuclc- 
ation of vortices. 

Motivated by both observations of coherent optical in- 
formation storage in atom clouds and the interesting dy- 
namics which slow and stopped light pulses can induce 
in BECs, the present paper theoretically examines new 
possibilities offered by stopping light in BECs. We first 
present a novel treatment of the switching process and 
establish that the switch-off and switch-on of the cou- 
pling field can occur without dissipating the coherent in- 
formation, provided the length scale of variations in the 
atomic wavefunctions are sufficiently large. We find that 
the switching of the coupling field can be done arbitrarily 
fast in the ultra-slow group velocity limit V„ -C c. This 
is consistent with previous work 0, 0, 0, , but goes 
further as it applies to pulses in both the weak and strong 
probe regimes [14| . We also present the first explicit cal- 
culation of how the behavior of the system differs when 
one switches the coupling field faster than the natural 



lifetime of the excited state T _1 (see Fig. ^a)) and find 
the information is successfully transferred even in this 
regime. Furthermore we will see the analysis presented 
here, being phrased in terms of the spatial characteristics 
of the atomic wavefunctions, is well suited to addressing 
the issue of storage for times long compared to the time 
scale for atomic dynamics. 

Investigation of this very issue of longer storage times 
constitute the central results of the paper. In this regime, 
there are important differences between optical informa- 
tion storage in BECs versus atom clouds above the con- 
densation temperature (thermal clouds), which have been 
used in experiments to date. During the storage time, the 
light fields are off and the atomic wavefunctions evolve 
due to kinetic energy, the external trapping potential, 
and atom- atom interactions. (We label all these external 
dynamics to distinguish them from couplings between the 
internal levels provided by the light fields.) To success- 
fully regenerate the probe pulse by then switching on 
the coupling field, there must be a non-zero value of the 
ensemble average (over all atoms) of the coherence be- 
tween |1) and 1 2). When revivals were attempted after 
longer storage times in thermal clouds (several millisec- 
onds in Q) the external dynamics had washed out the 
phase coherence between |1) and |2) (due to the atoms 
occupying a distribution of single particle energy levels) 
and no output pulses were observed. In zero temperature 
BECs the situation is completely different. There each 
atom evolves in an identical manner, described by a non- 
linear Schrodinger equation (Gross-Pitaevskii equation 
[To|). preserving the ensemble average of the coherence 
during the external dynamics. That is, even as the am- 
plitude and phase of the wavefunctions representing the 
BEC evolve, the relative phase of the two components 
continues to be well defined at all points in space. Thus, 
if we switch the coupling field back on, we expect that 
the evolved BEC wavefunctions will be written onto a 
revived probe field. In this way, the probe pulses can be 
processed by the BEC dynamics. 

We show several examples, relevant to current Rb-87 
BECs, of the interesting two component dynamics which 
can occur during the storage time. This atom possesses 
internal states with very small inelastic loss rates [l5|. 
and thus, long lifetimes We find that, depending 

upon the relative scattering lengths of the internal states 
involved and the relative intensity of the probe and cou- 
pling fields, one could observe the formation of interfer- 
ence fringes, gentle breathing motion, or the formation 
and motion of two-component (vector) solitons |17| . In 
particular we find that using the levels \F — 1, Mp = —1) 
and \F = 2, M F = +1) in Rb-87 could allow for long, ro- 
bust storage of information and controllable processing. 
In each case, we observe the amplitude and phase due 
to the dynamical evolution is written onto revived probe 
pulses. These pulses then propagate out of the BEC as 
slow light pulses, at which point they can be detected, 
leaving behind a BEC purely in its original state 

For practical applications of this technique to stor- 
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age and processing, one must understand in detail how, 
and with what fidelity, the information contained in the 
atomic coherence is transferred and output on the light 
fields. Thus, the last part of the paper is devoted to 
finding the exact relationship between the BEC wave- 
functions before the switch-on and the observed output 
probe pulses. We find equations linking the two in the 
ideal limit (without absorption or distortion). Then, us- 
ing our earlier treatment of the switching process, we 
identify several sources of imperfections and calculate the 
fidelity of the information transfer for Gaussian shaped 
wavefunctions of various amplitudes and lengths. In the 
weak probe limit, we find a simple relationship between 
the wavefunction in |2) and the output probe field. We 
find that optimizing the fidelity involves balancing con- 
siderations related to, on one hand, absorption of small 
length scale features in the BEC wavefunctions and, on 
the other hand, imperfect writing of wavefunctions which 
are too near the condensate edge. For stronger probes, 
we find a more complicated relationship, though we see 
one can still reconstruct the amplitude and phase of the 
wavefunction in |2) using only the output probe. In this 
regime, we find that the fidelity of the writing and out- 
put is nearly independent of the probe strength when 
the oscillator strengths for the two transitions involved 
(|1) «-» |3) and |2) <-> |3)) are equal. However, unequal os- 
cillator strengths lead to additional distortions and phase 
shifts, and therefore lower fidelity of information transfer, 
for stronger probes. 

We first, in Section^ introduce a formalism combin- 
ing Maxwell-Schrodinger and Gross-Pitaevskii equations 
|10| , which self-consistently describe both the atomic (in- 
ternal and external) dynamics as well as light field prop- 
agation. We have written a code which implements this 
formalism numerically. Section IIHI presents our novel 
analysis of the switching process, whereby the coupling 
field is rapidly turned off or on, and extends previous 
treatments of the fast switching regime, as described 
above. Section lTvl shows examples of the very rich variety 
of two-component dynamics which occur when one stops 
light pulses in BECs and waits for much times longer than 
the characteristic time scale for atomic dynamics. Sec- 
tion [V] contains our quantitative analysis of the fidelity 
with which the BEC wavefunctions are transferred onto 
the probe field. We conclude in Section IVII and antic- 
ipate how the method studied here could eventually be 
applied to transfer of information between distant BECs, 
and the generation of light with squeezed statistics. 



II. DESCRIPTION OF ULTRA-SLOW LIGHT 
IN BOSE-EINSTEIN CONDENSATES 

We first introduce our formalism to describe the sys- 
tem and review USL propagation within this formalism. 
Assume a Bose-condensed sample of alkali atoms, with 
each atom in the BEC containing three internal (elec- 
tronic) states in a A-configuration (Fig.QJa)). The states 



|1) and |2) are stable and the excited level |3), radiatively 
decays at T. In the alkalis which we consider, these in- 
ternal levels correspond to particular hyperfine and Zee- 
man sublevels for the valence electron. When the atoms 
are prepared in a particular Zeeman sub-level and the 
proper light polarizations and frequencies are used 0, |j| , 
this three level analysis is a good description in practice. 
All atoms are initially condensed in |1) and the entire 
BEC is illuminated with a coupling field, resonant with 
|2) «-* |3) transition and propagating in the +z direction 
(see Fig.^b)). A pulse of probe field, with temporal half- 
width To, resonant with the |1) <-» |3) transition, and also 
propagating in the +z direction, is then injected into the 
medium. The presence of the coupling field completely 
alters the optical properties of the atoms as seen by the 
probe. What would otherwise be an opaque medium 
(typical optical densities of BECs are ~ 400) is rendered 
transparent via electromagnetically induced transparency 
(EIT) Q, and the light pulse propagates at ultra-slow 
group velocities (~ 10 m/s)@. As this occurs, all but a 
tiny fraction of the probe energy is temporarily put into 
the coupling field, leading to a compression of the probe 
pulse to a length smaller than the atomic medium itself 



To describe the system theoretically, we represent 
the probe and coupling electric fields with their Rabi- 
frequencies £l p ( c ) = — d 13 ( 2 3) 1 E p ( c )/7l, where E p ( c j are 
slowly varying envelopes of the electric fields (both of 
which can be time and space-dependent) and d 13 ( 2 3) are 
the electric dipole moments of the transitions. The BEC 
is described with a two-component spinor wavefunction 
("01 j "02 ) T representing the mean field of the atomic field 
operator for states |1) and |2). We ignore quantum fluc- 
tuations of these quantities, which is valid when the tem- 
perature is substantially below the BEC transition tem- 
perature fl9[ . The excited level |3) can be adiabatically 
eliminated [2(J when the variations of the light fields' en- 
velopes are slow compared to the excited state lifetime 
r _1 (which is 16 ns in sodium). The procedure is out- 
lined in the Appendix. The functions ?pi , if)% evolve via 
two coupled Gross-Pitaevskii (GP) equations jlG]. For 
the present paper, we will only consider dynamics in the 
z dimension, giving the BEC some cross-sectional area 
A in the transverse dimensions over which all dynamical 
quantities are assumed to be homogeneous. This model 
is sufficient to demonstrate the essential effects here. We 
have considered effects due to the transverse dimensions, 
but a full exploration of these issue is beyond our present 
scope. The GP equations are [iSl l2ll]: 
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where m is the mass of the atoms and we will consider a 
harmonic external trapping potential Vi(z) = ^muj z 2 z 2 . 
The potential V% can in general differ from V\. For ex- 
ample, in a magnetic trap the states |1) and |2) can have 
different magnetic dipole moments. However, in the ex- 
amples we consider the potentials are equal to a good 
approximation, and so V2(z) — V\{z) is assumed in the 
calculations. Atom-atom interactions are characterized 
by the Utj — 4irN c h 2 aij /mA, where N c is the total num- 
ber of condensate atoms and are the s-wave scattering 
lengths for binary collisions between atoms in internal 
states \i) and \j). The last pair of terms in each equa- 
tion represent the coupling via the light fields and give 
rise to both coherent exchange between |1), |2) as well as 
absorption into |3). In our model, atoms which populate 
1 3) and then spontaneously emit are assumed to be lost 
from the condensate (Fig. 0a)), which is why the light 
coupling terms are non-Hcrmitian. 

The light fields' propagation can be described by 
Maxwell's equations. Assuming slowly varying envelopes 
in time and space (compared to optical frequencies and 
wavelengths, respectively) and with the polarization den- 
sities written in terms of the BEC wavefunctions these 
equations are [fllElj: 



d_ ld_\ 

dz + cdt) p 

d 1 d \ _ N c f 2 3<7o /r> 



.^p(fi p |^| 2 +fi c ^V 2 ), 



dz cdt c 2 A 



(Jlclifel +^iV>2 12) 



where /i3,/23 are the dimensionless oscillator strengths 
of the transitions and <jq is the resonant cross section 
(1.65 x lCT 9 cm 2 and 2.91 x l(T 9 cm 2 for the D 1 lines of 
sodium and Rb-87, respectively). These equations ignore 
quantum fluctuations of the light fields, just as the GP 
equations ignore quantum fluctuations of the atomic 
fields. Note that in the absence of any coherence ipi^ 
these equations predict absorption of each field with the 
usual two-level atom resonant absorption coefficient (the 
atom density times the single atom cross-section). 

For much of our analysis, we will solve Eqs. 0-0 
self consistently with a numerical code, described in 
iEHH. In 0, this was successfully applied to pre- 
dict and analyze the experimental observation of the nu- 
cleation of solitons and vortices via the light roadblock. 
When the light fields are off, the two states do not 



exchange population amplitude and Eqs. reduce to 
coupled Gross-Pitaevskii equations for a two-component 
condensate. By contrast, when they are on and the probe 
pulse length To (and inverse Rabi frequencies fi" 1 , Qj 1 ) 
are much faster than the time scale for atomic dynamics 
(typically ~ ms), the internal couplings induced by the 
light field couplings will dominate the external dynamics. 

For our our initial conditions, we consider N c ~ 10 6 
atoms are initially in |1) in the condensed ground state. 

We determine the wavefunction of this state ip{ (z ) 
numerically by propagating in imaginary time [23j, 
though the Thomas-Fermi approximation [2J| provides 
a good analytic approximation to ip\ (z). In the cases 
presented in this paper, we choose the transverse area 
A so that the central density and chemical potential p, = 



(G), 



are in accordance with their values in a trap 



with transverse frequencies io,, 



3.8 u) z , as in pre- 
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vious experiments [!j ■ In Fig. |3 we consider an example 
with N c — 1.2 x 10 6 sodium atoms (with an = 2.75 nm 
|2flj). cv z = (2tt)21 Hz and A = tt(8.3 fim) 2 . The ground 
state density profile N c \ip^ (z)\ 2 is indicated with the 
dotted curve in Fig.[5{b). In this case the chemical po- 
tential fi = (2?r) 1.2 kHz. 

With our initial ground state tp[ (z) determined, con- 
sider that we initially (t = —oo) input a constant cou- 
pling field with a Rabi frequency Q c q and then inject a 
Gaussian shaped probe pulse at z- la (see Fig. ^b)) with 
a temporal half- width to and a peak Rabi frequency Sl p o- 
We define our times such that t = corresponds to the 
time the peak of the pulse is input. The dotted and 
dashed curves in Fig. Ufa) show, respectively, a constant 
coupling input Sl c0 = (2tt) 8 MHz and a weaker input 
probe pulse with peak amplitude Q p o — (27r)3.5 MHz and 
width To = 1.5 /is. 

Solving Eqs. 0-© reveals that, when the conditions 
necessary for EIT hold, the pulse will compress upon en- 
tering the BEC and propagate with a slow group velocity. 
As it does this, it transfers the atoms into superpositions 
of |1) and 1 2) such that plf : 



ip2(z,t) 



Q p (z,t) 

'n c ( z ,t) 



ipi{z,t), 



(3) 



which is a generalization to BECs of the single atom dark 
state |26| . Thus, the probe field Sl p imprints its (time 
and space dependent) phase and intensity pattern in the 
BEC wavefunctions as it propagates. When © is ex- 
actly satisfied, it is easily seen that the two light field 
coupling terms in each of Eqs. cancel, meaning that 
|1) and 1 2) are completely decoupled from the excited 
state 1 3). However, time dependence of f2 p causes small 
deviations from Q to occur, giving rise to some light- 
atom interaction which is, in fact, the origin of the slow 
light propagation. If one assumes the weak probe limit 
(f2 p o <IC fi c o) and disregards terms of order (f2 p o/^co) 2 
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FIG. 2: Coherent storage of a light pulse via fast 
switching. Results of numerical integration of Eqs. Q- 
@ showing both a slow and stopped light experiment, (a) 
The dashed curve shows the normalized input probe pulse 
intensity |f2 p (zi n ,i)| 2 /fipo ( a r o = 1.5 [ts 1/e half-width, 
Q p o = (27r) 3.5 MHz pulse), while the dotted curve shows 
|fic(2m,t)| 2 /fico f° r a constant input coupling field Q c o = 
(2ir) 8.0 MHz. The thick solid curve shows the delayed out- 
put probe pulse (at z out , see Fig. [Jb)). The time that the 
tail of the input has vanished and the rise of the output has 
not begun (marked by the arrow) corresponds to the time 
the pulse is completely compressed inside the BEC. In the 
output coupling intensity (thin solid curve) we see adiaba- 
tons. The values cr = 1.65 x 10~ 9 cm 2 and /13 = 1/2 
and /23 = 1/3 have been used, (b) The atomic density 
in the two states, N c \il:i\ 2 /A (solid curve) and N c \i/j2\ 2 /A 
(dashed curve), at t = 3.7 (ms (indicated by the arrow in 
(a)). The dotted curve shows the orig inal density N c \i/)[ G) \ 2 /A 
before the probe is input. The arrow here and in (c)-(d) 
marks the position z c discussed later in the text, (c) Spatial 
profiles of the probe r2 p | 2 /f2 2 (solid curves) and coupling 
|f2 c | 2 /f2 2 (dashed curves) field intensities at various times 
while the coupling field input is switched off. The switch-off 
is an error function profile with a width t s — 0.1 [is, cen- 
tered at tog = 3.7 /MS. Successively thinner curves refer to 
t = 3.52, 3.65, 3.68, and 3.75 /is. The dotted curve indicates 
the original condensate density (arb. units), (d) Spatial pro- 
files of the normalized dark |S7.d | 2 /f2 2 (dashed curves, scale 
on the left) and absorbed |f2>i| 2 /f2 2 o (solid curves, scale on 
right) field intensities, at the same times as in (c). (e) The 
input and output probe and coupling intensities, with the 
same conventions as (a), in a stopped light simulation. In 
this case, the coupling field is switched off at t a g = 3.7 pis and 
then back on at ton = 15.3 fis. 

the group velocity of the probe pulse is 0, : 



and so is proportional the coupling intensity and inversely 



proportional to the atomic density N c \ip[ (z)\ 2 /A. The 
half-width length of the pulse in the medium is L p — 
T oV g , which is a reduction from it's free space value by 
a factor V g /c, while its peak amplitude does not change. 
Thus only a tiny fraction of the input probe pulse energy 
is in the probe field while it is propagating in the BEC. 
Most of the remaining energy coherently flows into the 
coupling field and exits the BEC, as can be seen in the 
small hump in the output coupling intensity (thin solid 
curve in Fig. |2Ja)) during the probe input. This has 
been dubbed an adiabaton 27]. The magnitude of the 
adiabaton is determined by the intensity of the probe 
pulse. 

In the weak probe limit ipi(z,t) never significantly 
deviates from the ground-state wavefunction ipi(z,t) ~ 
V4 (z) and the coupling field is nearly unaffected by the 
propagation Q c (z,t) « Q c o (both of these results hold to 
O(£l 2 Q /il 2 )). I n t ms casc © shows tp2 follows the probe 
field tl p as the pulse propagates. The arrow in Fig. |2fa) 
indicates a time where the probe has been completely in- 
put and has not yet begun to output. During this time, 
the probe is fully compressed in the BEC and Fig. |2f b) 
shows the atomic densities in |1) and |2) at this time. The 
spatial region with a non-zero density in \2) corresponds 
to the region occupied by the probe pulse (in accordance 
with (J3J). Equation J2| applies to the phases as well as 
the amplitudes, however the phases in this example are 
homogenous and not plotted. 

Once the pulse has propagated through the BEC, it 
begins to exit the +z side. The energy coherently flows 
back from the coupling to probe field, and we see the 
output probe pulse (thick solid curve). Correspondingly, 
we see a dip in the coupling output at this time. In 
the experiments the delay between the input and output 
probe pulses seen in Fig. Ufa) is measured with a PMT 
(Fig.[IJb)). This delay and the length of the atomic cloud 
is used in to calculate the group velocity. The group 
velocity at the center of the BEC in the case plotted is 
y fl (0)=6m/s. 

Note that the output pulse plotted in Fig. Eta) at z ou t 
is slightly attenuated. This reduction in transmission is 
due to the EIT bandwidth. The degree to which the 
adiabatic requirement ro T/Sl 2 Q is not satisfied, will 
determine the deviation the wavefunctions from the dark 
state ©, which leads to absorption into |3) and subse- 
quent spontaneous emission. Quantitatively, this reduces 
the probe transmission (the time integrated output en- 
ergy relative to the input energy) to |21| : 

(5) 

where D{z) = (N c f 13 (r a /A) f* n dz'\ip[ G) (z')\ 2 is the op- 
tical density. In the Fig. example, D(z ou t) = 390. 
The peak intensity of the pulse is reduced by a factor 
T 2 while the temporal width is increased by T _1 (this 
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spreading can be seen in Fig. Ufa)). The appearance 
of the large optical density D(z out ) in JSJ represents 
the cumulative effect of the pulse seeing a large num- 
ber of atoms as it passes through the BEC. To prevent 
severe attenuation and spreading, we see we must use 

to > Tq" 11 "' = 2y/D(z ut) r/^cOi which is 1.0 fis in our 
example. 



III. FAST SWITCHING AND STORAGE OF 
COHERENT OPTICAL INFORMATION 

We now turn our attention to the question of stop- 
ping, storing and reviving probe pulses. We will show 
here that once the probe is contained in the BEC, the 
coupling field can be switched off and on faster than the 
EIT adiabatic time scale without causing absorptions or 
dissipation of the information 0, . While we find no 
requirement on the time scale for the switching, we will 
obtain criteria on the length scales of "01 , 02 which must 
be maintained to avoid absorption events. While pre- 
vious work 0, IT2I IT^ | had addressed the fast switching 
case in the weak-probe limit, here we obtain results that 
are valid even when f2 p o ~ ^co UM- As we will see, for 
shorter storage times where the external atomic dynam- 
ics play no role, these length scale requirements can be 
related to the adiabatic requirements on the input probe 
pulse width To. The advantage of the current analysis is 
that it is easily applied to analyzing probe revivals after 
the atomic dynamics have completely altered the wave- 
functions ip\,ip2- In such a case it is these altered wave- 
functions, rather than the input pulse width To, which is 
relevant. 



Analyzing fast switching in the dark/absorbed 
basis 



Consider that we switch the coupling field input off at 
some time t g with a fast time scale t s -C tq. Fig. Etc) 
plots the probe and coupling intensities as a function of z 
at various times during a switch-off with t s =0.1 /is. We 
see the probe intensity smoothly ramps down with the 
coupling field such that their ratio remains everywhere 
constant in time. Remarkably for reasons we discuss 
below, the wavefunctions "01 > V'2 are completely unaffected 
by this switching process. Motivated by this, we will, in 
the following, assume that V'liV^ do not vary in time 
during this fast switching period, and later check this 
assumption. 

To understand this behavior, it is useful to go into a 
dark/absorbed basis for the light fields, similar to that 
used in [H, by defining 
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-r 2 n 

01 4>2 



where ipo = V^0iF + IV^I 2 - From JBJ) one sees that when 
the condensate is in the dark state 0, &a = 0. Using the 
notation ipi = |^|e Ic ^ and transforming the propagation 
equations © according to iJBJ, one gets: 
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and we have ignored the vacuum propagation terms 
~ 1/c in Eq. They are unimportant so long as the 
fastest time scale in the problem is slow compared to 
time it takes a photon travelling at c to cross the conden- 
sate (about 100 fim/c ~ lps). Note that a A represents 
the usual absorption coefficient weighted according to the 
atomic density in each of |1) and |2). The terms 0£n A , oti 
arise from spatial variations in the wavefunctions, which 
make the transformation © space-dependent. The term 
ct\i represents an additional effect present when the light- 
atom coupling coefficient differs on the two transitions 
(/i3 /23) and is discussed in detail in Section IV CI 

Consider for the moment a case with /13 = /23 (im- 
plying ai2 = 0) and assume a region in which ipi,ip2 are 
homogenous (so aN A = a>i = 0). Equation © shows 
that the absorbed field £l A attenuates with a length scale 
a~ A x , the same as that for resonant light in a two- level 
atomic medium, and less than 1 /xm at the cloud center 
for the parameters here. One would then get £l A — > 
after propagating several of these lengths. Conversely, 
shows the dark light field £1 d experiences no interac- 
tion with the BEC and propagates without attenuation 
or delay. 

However, spatial dependence in 0i,"02 gives rise to 
otNA,cti 7^ in 0-©, introducing some coupling be- 
tween Qd and Q A , with the degree of coupling governed 
by the spatial derivatives dipi/dz, dip2/dz. A simple and 
relevant example to consider is the case of a weak ultra- 
slow probe pulse input and contained in a BEC, as dis- 
cussed above (see Fig. IHb)). The wavefunction -02 has 
a homogenous phase and an amplitude which follows the 
pulse shape according to ©, meaning that olna scales 
as the inverse of the pulse's spatial length L~ l , 

It is important to note that this coupling is determined 
by spatial variations in the relative amplitude "02 /ipi and 
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does not get any contribution from variations in the to- 
tal atomic density ipQ. To see this we note that if we 
can write the wavefunctions as 0i(z) = Ciipo(z), ip2(z) = 
C2tpo(z), where c\, C2 are constants independent of z, then 
a at a and a\ as defined in © vanish. 

When variations in the relative amplitude are present, 
examination of (jSJ) reveals that when the damping is 
much stronger than the coupling (|ajvA|,c*z <C a a), one 
can ignore the spatial derivative term in analogy to an 
adiabatic elimination procedure. When we do this, © 
can be approximated by: 



una 

OiA 



(10) 



Strictly speaking, one can only apply this procedure in 
the region where £Ia has propagated more than one ab- 
sorption length (that is for z > z c where z c is defined 
by J z c dz'aA(z') = 1). However, z c is in practice only 
a short distance into the BEC (marked with arrows in 
Figs.|3b)-(d)). When the probe has already been input, 
as in Fig. 2(b), ip2 and therefore £Ia (see ©) are already 
trivially zero for z < z c so 111 Oft holds everywhere. We 
next plug Eq. (fTU|) into Eq. (|7|), giving: 
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These last two equations rely only on assumptions about 
the spatial derivatives of ip\ , "02 and not on the time scale 
of the switch-off t s (though we have assumed t s 3> in 
adiabatically eliminating |3) from our original Eqs. 
0). 

These results allow us to conclude two important 
things which hold whenever |ajvA|,c*2 a a and the 
probe has been completely input. First, the coefficients 
governing the propagation of flp in 1|11|) are extremely 
small. The length scales IowaI" 1 , oi[2 are already gen- 
erally comparable to the total BEC size and the length 
scales for changes in 91 d, given by (fTTJ) . scale as these 
terms multiplied by the large ratio ola/chna- (The ai 
term in (|llfl can lead to additional phase shifts, which 
we discuss in Section IV Cfl . Therefore the dark field 
9lo propagates with very little attenuation. As we have 
noted (and later justify and discuss in more detail) the 
wavefunctions (and therefore the propagation constant 
in brackets in (|llfl ) are virtually unchanged during the 
switch-off. Under this condition, changes initiated in 91d 
at the entering edge z;„ quickly propagate across the en- 
tire BEC. To apply this observation to the switch-off, we 
note when the pulse is contained so 02 = at the entering 
edge, © shows ftp — 9l c there. Switching off the cou- 
pling field at Zi n then amounts to switching off 91 p> at z- ln 
and (|llfl shows that this switch-off propagates through 
the entire BEC with little attenuation or delay. Second, 
from Q10[) we see that as 9lrj is reduced to zero ft a is 
reduced such that the ratio 91 a /^d remains constant in 
time. 



A numerical simulation corroborating this behavior is 
plotted in Fig. |2Id). The dark field intensity \9lo\ 2 is 
seen to switch-off everywhere as the coupling field in- 
put is reduced to zero over a r s = 0.1 /is timescale, con- 
firming that the changes in the input propagate across 
the BEC quickly and with little attenutation. The ra- 
tio |r2^| 2 /|f2£>| 2 is everywhere much smaller than unity 
and constant in time. The only exception to this is in a 
small region z < z c at the cloud entrance, where 91a has 
not yet been fully damped. The plot of \91a\ 2 during the 
switch-off indeed demonstrates how (|10[) is a generally 
good approximation, but breaks down in this region. In 
the case plotted (and any case where the probe is fully 
contained) the wavefunction 02 is so negligible in this re- 
gion that 91a is rather small and unimportant (note the 
scale on the right hand side of the plot). 

Translating this back into the 9l p , 9l c basis, we note 
from (|SJ| that keeping 91a ~ means the probe must (at 
all z) constantly adjust to the coupling field via: 



a, 



(12) 



Thus we see that 9l p smoothly ramps down with 9l c even 
if 9l c is ramped down quickly, as seen in Fig.[2Ic). 

Using our results for the light fields, we can now see 
why the wavefunctions ipi , ip2 do not change during the 
switching. Physically, the probe is in fact adjusting to 
maintain the dark state H12|) . and in doing so induces 
some transitions between |1) and |2). However, only a 
fraction V g /c ~ 10~ 7 of the input energy is in the probe 
while it is contained in the medium, and so the probe is 
completely depleted before any significant change occurs 
in 01,02 0- In fact, the energy content of the probe 
field right before the switch-off is less than 1/ 100th of 
a free-space photon in the case here. Thus, 9l p is com- 
pleted depleted after only a fraction of one |1) — > |2) 
transition. Note that H12fl is equivalent to ©■ However, 
writing it in this way emphasizes that, during the switch- 
ing, the probe is being driven by a reservoir consisting of 
the coupling field and atoms and adjusts to establish the 
dark state. This is contrast to the situation during the 
probe input, when many photons from both fields are be- 
ing input at a specific amplitude ratio, forcing the atomic 
fields to adjust to establish the appropriate dark state. 
Plugging in our results for the light fields l|10l) - (|llfl into 
our Eqs. we can calculate the changes that occur in 
■01 and i\>2 during the switch-off. Doing this, we find rela- 
tive changes in 0i, ip2 are both smaller than t s /tq, which 
can be made arbitrarily small for fast switching t s < To. 
The little change which does occur is due not to any pro- 
cess associated with the switching itself but is due to the 
small amount of propagation during the switch-off. It 
is therefore safe to assume (as we saw numerically) that 
0i , 02 are constant in time during the switch-off. One 
can also show with this analysis that the ratio of co- 
herent exchange events (from |1) to |2) or vice- versa) to 
absorptive events (transitions to |3) followed by sponta- 
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neous emission) is |f2_D/f2;i| implying that the switch-off 
occurs primarily via coherent exchanges. 

In the original stopped light experiment 3, the cou- 
pling field was then switched back on to f2 c o, after some 
controllable storage time T s t , at a time t on — i ff + T s t . In 
that experiment, revivals were observed for storage times 
r st too short for significant external atomic dynamics to 
occur, so the state of the atoms was virtually identical 
at t on and t a g. Then the analysis of the switch-on is 
identical to the switch-off as it is just the same coher- 
ent process in reverse. The probe is then restored to the 
same intensity and phase profile as before the switch-off. 
An example of such a case is plotted in Fig. Efe). The 
output revived pulse then looks exactly like the normal 
USL pulse (compare the output in Figs. [Ha), (e)), as was 
the case in the experiment. 

We have established the requirement Ioat^I -C a a 
is necessary and sufficient for coherent switching to oc- 
cur. Under what conditions is this satisfied? When the 
switch-off occurs while the probe is compressed, it is al- 
ways satisfied because of bandwidth considerations men- 
tioned above (see Eq. |J5J). Specifically, the input pulse 



must satisfy tq > Tq 



(min) 



However, this leads to a pulse 



width in the medium of L p > 2^ D(z out ) a^ 1 , implying 
that \o!na\ "C a a is satisfied (as D(z out ) » 1). There- 
fore, any pulse which can successfully propagate to the 
cloud center, can be abruptly stopped and coherently 
depleted by a rapid switch-off of the coupling field. Of 
course, if the switch-on is then done before significant 
atomic dynamics, the same reasoning applies then. In 
this case, our requirements on the spatial derivatives are 
already encapsulated by the adiabatic requirements on 
the probe pulse. It is when the external atomic dynam- 
ics during the storage significantly change ipi,tp2, that 
the analysis of the switch-on becomes more complicated. 
This is the central purpose of the Sections IIVIVI below. 

As hinted above, when the pulse is not yet completely 
input, the switching can cause absorptions and dissipa- 
tion of the information. In this case tp2 is significant in 
the region z < z c and so Qa is significant for several ab- 
sorption lengths into the BEC. Physically, the coupling 
field sees atoms in a superposition of |1) and |2) imme- 
diately upon entering the condensate rather than only 
|1) atoms. Numerical simulations of this situation con- 
firm this. During both the switch-on and switch-off, a 
significant number of atoms are lost from both conden- 
sate components, primarily concentrated in the z < z c 
region, and the revived probe pulse is significantly atten- 
uated relative to the pulse before the storage. 

Incidentally, this explains the apparent asymmetry be- 
tween the probe and coupling in a stopped light experi- 
ment. Whenever both fields are being input, the temporal 
variations in both fields must be slow compared with the 
EIT adiabatic time scale, as the atomic fields must adjust 
their amplitudes to prevent absorptions. However, when 
one of the fields is no longer being input (like a contained 
probe pulse), the other field's input can be quickly varied 
in time. 
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FIG. 3: Switching faster than the natural linewidth. 

Solid and dashed lines show, respectively, the normalized 
probe and coupling amplitudes at a point in the center of 
the BEC Q p (z = -12 /im,t)/(] p0 , M c (z = -12 fj,m,t)/fl c0 
when the switching is faster than the natural linewidth r s = 
2 ns <C r _1 (but otherwise the same parameters as Fig.|2fe)). 
The two insets magnify the regions near the t a g and t on . 



B. Switching faster than F 1 

When the switching is done quickly compared to the 
natural lifetime of |3) (t s < T _1 ), the assumptions 
needed for adiabatic elimination of ip3 are no longer valid. 
The analysis of the switching then becomes more compli- 
cated but, remarkably, we find the quality of the storage 
is not reduced in this regime. 

To see this, we performed numerical simulations of 
switching in this regime by solving Eqs. (|A.1(I - |IA.4I) 
(which arc the analogues of Q-lJlJl without the adiabatic 
elimination of -03) for the same parameters as Fig.[2]but 
varying r s from 80 ns to 0.1 ns 29|. An example with 
T,,. = 2 ns is shown in Fig. We plot the probe and 
coupling amplitudes versus time at a point near the cen- 
ter of the BEC, and zoom in on the regions near the 
switch-on and -off. Note that that while the coupling 
field smoothly varies to it's new value in each case, the 
probe field slightly overshoots the value given by i|12|) be- 
fore returning to it. In all cases with r s < the probe 
amplitude \Q p \ experienced damped oscillations before 
reaching its final value. The frequency of the oscillations 
was determined by t~ and the damping rate was F. 

During these oscillations, the absorbing field £Ia is in 
fact quite significant. However, the maximum value of 
| ^3 1 and the total atomic loss due to spontaneous emis- 
sions events was completely independent of t s over the 
range we explored. To understand this, we note that 
even though Qa is larger, the frequency components of 
the quickly varying coupling field are spread over a range 
much wider than T. Thus a very small fraction of the 
energy in the field f2 c is near resonant with the atomic 
transition and it passes through the BEC unattenuated. 
Analyzing the problem in Fourier space one can see that 
the increase in £Ia for smaller r s is precisely offset by this 
frequency width effect, leading us to the conclusion that 
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spontaneous emission loss is independent of t s . 

Unlike the case of switching slowly compared with , 
the adjustment of the probe is not a coherent process. 
Rather, spontaneous emission damps the system until the 
dark state (|12[1 is reached. However, because we are in 
the regime < c the amount of energy in the probe is 
much less than one photon, and thus still has virtually 
no impact on tb-\ , ife. This will not necessarily hold when 
V g ~ c and in [l3| it was predicted that there are, in fact, 
adiabatic requirements for t s in this regime. However, in 
all cases of interest here this inequality is well satisfied 
and there is no restriction on r s . 



IV. BEC DYNAMICS AND PROCESSING 
OPTICAL INFORMATION 

We now turn to the question of the atomic dynam- 
ics during the storage time. These dynamics depend 
strongly on the relative scattering lengths of the states 
used, and the probe to coupling intensity ratio. We will 
see that the these dynamics are written onto revived 
probe pulses by switching the coupling field back on. 
The length scale requirement on variations in the rela- 
tive amplitude ip2/^i derived above (|ajvA| a a) plays 
a central role in the fidelity with which these dynamics 
are written. We will demonstrate with several examples, 
relevant to current Rb-87 experimental parameters. 



A. Writing from wavefunctions onto light fields 

■ ■ 

In [4j the revived pulses were seen to be attenuated 
for longer storage times with a time constant of ~ 1 ms. 
The relative phase between the |1) and |2) (or coher- 
ence), averaged over all the atoms, washed out because 
the atoms, although at a cold temperature of 0.9 fiK, 
were above the Bose-condensation temperature and so 
occupied a distribution of energy levels. By contrast, a 
zero temperature two-component BEC will maintain a 
well defined phase (the phase between ipi and 1P2) at all 
z, even as ?Ai,'02 evolve. They evolve according to the 
two-component coupled GP equations Q with the light 
fields set to zero (f2 p = fi c = 0), and the initial condi- 
tions are determined by the superposition created by the 
input pulse at i ff (see Pig. E[b)). After such an evolu- 
tion, we can then switch the coupling field back on and 
the probe will be revived according to (|12|) but with the 
new wavefunctions ipi,ipz, evaluated at t on . 

Here we examine a wide variety of different two- 
component BEC dynamics which can occur after a pulse 
has been stopped. Because of the initial spatial structure 
of ipi,ip2 created with the USL pulses, we see some novel 
effects in the ensuing dynamics. 

Upon the switch-on, in many cases the probe will be 
revived via (|12[1 . We label the spatial profile of the re- 
vived probe by Q p (2) = O p (z, t on + t s ). When we are 



in the weak probe limit (|f2 p | <C |O c | or, stated in terms 
of the wavefunctions, \ip2\ <C |V>i|), becomes 

n(-) ( , ) = _V^) Oc0 . (13) 

allowing us to make a correspondence between the re- 
vived probe and ip2 if f2 c o and ip\ (z) are known. The 
revived pulse will then propagate out of the BEC to the 
PMT at z out . In the absence of any attenuation or distor- 
tion during the propagation out, the spatial features of 
the revived probe get translated into temporal features. 
Thus if we observe the output 17 p out '' (t) = f2 p (z out , t) then 
we would deduce that revived pulse was 

r» p rov ) (z) = r» p out ) (r(z out ) - t(z) + t on ) (14) 

where r(z out ) - r(z) = (T/n 2 c0 )(D{z out ) - £>(*)) is the 
time it takes the probe pulse to travel from some point 
z to O ut- Note that we are relying here on the fact that 
we can switch the coupling field onto its full value of Q c q 
with a time scale fast compared to the pulse delay times 
T~(z ou t) — T ( z )- A slow ramp up would lead to a more 
complicated relationship between f2 p r (z) and f2 p out ' > (t). 

Combining and l|14|l shows how the phase and 
amplitude information that was contained in ip2 at the 
time of the switch-on is transferred to the output probe 
£lp {t). This transfer will be imperfect for three rea- 
sons. First, when sufficiently small spatial features are in 
■02 j then ckata is comparable to a a, giving rise to a signifi- 
cant Qa- The resulting absorptions will cause deviations 
of Cl p ley ' from our expectation l|13|) . Second, we men- 
tioned how spatial features are translated into temporal 
features on the output probe. During the output, fast 
time features on the output probe will be attenuated via 
the bandwidth effect discussed in J5J , affecting the accu- 
racy of the correspondence 114|) . Finally, stronger output 
probes, which will occur when \ip2 \ ~ IV'il at the time of 
the switch-on, make both the writing at the switch-on 
and the subsequent propagation out more complicated 
and less reliable. Our following examples will demon- 
strate these considerations. 



B. Formation and writing of interference fringes 

Interesting dynamics occur in when the two internal 
states are trapped equally (V2 = Vi) and the scatter- 
ing lengths ai2,an are slightly different. We consider 
a case with N c = 1.0 x 10 6 Rb-87 atoms and choose 
|1) = |55 1/2 ,F = 2,M F = +1), \2) - \5S 1/2 ,F = 
1,M F = -1), and 1 3) = \5P 1/2 ,F = 2,M F = 0). The 
two lower states |1) and |2) are magnetically trapped 
with nearly identical magnetic moments, and we con- 
sider a trap with lv z = (2tt) 21 Hz and assume a tran- 
vserse area A — 7r(5 /Ltm) 2 . These two states have an 
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FIG. 4: Formation of interference fringes in Rb-87 and 
writing fringes to probe light field, (a) The initial wave- 
functions created by a r = 2.0 [is, Q p o = (2n) 2 MHz probe 
pulse. The coupling field is initially on at Q c o = (2tv) 8 MHz 
then switched off at t Q g — 7.5 /is. The cloud parameters 
and atomic levels are described in the text. The thick solid 
curve shows the density N c \ip2(z)\ 2 /A (multiplied by 10 for 
clarity) while the thick dashed curve shows the phase i>i{z). 
The thin curve shows the density N c \ipi(z)\ 2 /A. The phase 
4>\ is approximately homogenous and not shown, (b) The 
effective potential V2-off//i l|15fl seen by the condensate in 
|2). (c) The resulting evolution of the density N c \ip2(z,t)\ 2 / A 
(white represents higher density) as the wave packet is pushed 
down the both sides of potential hill and interferes with it- 
self upon reaching the hard wall in V^-eff at the condensate 
edge. The vertical grey lines show the condensate bound- 
ary, and the horizontal grey lines indicate the times plotted 
in (d). (d) Plots of the density and phase of the wavefunc- 
tions at the indicated times (with the same conventions as 
(a)) at t = 55 ms, 101 ms and 291 ms. (e) Spatial profiles 
of the probe and coupling field directly after a fast switch- 
on after a storage time r st = t on — t a g corresponding to the 
times in (d) . The thick (thin) solid curves show the intensities 
\n ( p cv> {z)/n p0 \ 2 {\n c (z)/Q c0 \ 2 ) and the dashed curve shows 
the phase (f> p lev ^ (<f> c is nearly homogenous and not plotted), 
(f ) The time profile of the probe output intensity (solid curve) 
and phase (dots) at the output z out . 



anomalously small inelastic collisional loss rate [15J and 
so have been successfully used to study interacting two- 
component condensates for hundreds of milliseconds . 
The elastic scattering len gths are an = 5.36 nm, a\ 2 = 
1.024an, a 22 = 1.057an ^3. 

Figure 0f a) shows the initial wavefunctions tpi , ip 2 af- 
ter the input and stopping of a probe pulse. We see a 
Gaussian shaped density profile in |2) reflecting the in- 
put probe's Gaussian shape. The subsequent dynamics 



are governed by Eqs. Q with il p = Sl c = 0. In this 
case, a weak probe (Qp Q = S7^ /16) was input. As a re- 
sult, -01 (z, t) w tp^' (z) is nearly constant in time and the 
evolution of ip 2 is governed by essentially linear dynam- 
ics, with a potential determined by the magnetic trap 
and interactions with |1) atoms: 

V 2 -efi(z) = V 2 (z) + U 12 \iP[ G) (z)\ 2 - ii. (15) 

Figure E^b) shows this potential in this case. The hill 
in the middle has a height n[(ai 2 — an)/an] and arises 
because atoms feel a stronger repulsion from the conden- 
sate in |1) when they are in |2). Figures Elc-d) show 
the subsequent dynamics. One sees that the |2) conden- 
sate is pushed down both sides of the potential hill and 
spreads. However, once it reaches the border of the BEC, 
it sees sharp walls from the trap potential. Even for the 
fairly moderate scattering length difference here, there is 
sufficient momentum acquired in the descent down the 
hill to cause a reflection and formation of interference 
fringes near the walls. The wavelength of the fringes is 
determined by this momentum. 

What happens if one switches the coupling field back 
on after these dynamics? Figures 0Je) shows the revived 

probe pulses f2p rev ^ (z) upon switch-ons at the times cor- 
responding to Fig. 0fd) . One sees a remarkable transfer 
of the sharp density and phase features of the |2) con- 
densate onto the probe field, according to (|T3)l . Because 
we are in the weak probe regime, the coupling field in- 
tensity only very slightly deviates from it's input value 
Cl C Q. Figures^f) then shows the output probe ^^^(t). 
The sharp interference fringes are able to propagate out, 
though there is some attenuation and washing out of the 
features. Note that the features at more positive z prop- 
agate out first, leading to a mirror image- like relationship 
between the spatial Fig.^Je) and temporal Fig.^f) pat- 
terns as predicted by (|14f> . 

In this particular example, we successfully output 
many small features from the i\) 2 to the output probe 
field. Here the fringes are ~ 3 /im, which is still larger 
than the absorption length ~ 0.5 /im, but not substan- 
tially so, leading to some small amount of dissipation 
during the switch-on and output. This gives us a sense 
of the "information capacity". The number of absorp- 
tion lengths, or optical density, which in this case is 
-D(zout) ~ 300, ultimately limits the number of features 
(for a given desired fidelity) which could be successfully 
written and output. Note also that in the r st = 101 ms 
case, some of the ip 2 amplitude occupies the entering re- 
gion (z < z c ) leading to additional imperfections in the 
writing process in this region. 

C. Breathing behavior and long storage 

The dynamics in the previous example are quite dra- 
matic, but are not particularly conducive to preserving 





FIG. 5: Breathing behavior and long storage In this 
case the states chosen for |1) and \2) are reversed, changing 
the sign of the curvature in the effective potential 1151 in (a). 
All plots are the plotted with the same conventions as Fig. [I] 
except we additionally plot thin dotted curves in (c) showing 
the phase (f>i{z) and in (d) showing the phase </> c (z). 



or controllably processing the information in the BEC. 
For this, it would be preferable to switch the roles of 
|1) and |2). Such a case is shown in Fig. EJ Because 
"12 < a n in this case, the potential hill is turned into 
a trough (see Fig. OJa)). The effective potential in the 
region of the |1) condensate, in the Thomas-Fermi limit, 
becomes V^-eS — V\[(an — <2i2)/oui] and so it harmonic, 
with a much smaller oscillator frequency than the mag- 
netic trap. The evolution can be easily calculated by 
decomposing the wavefunction ip2 into a basis of the har- 
monic oscillator states of this potential. In the example 
here, there is significant occupation of the first several 
oscillator levels, and so one sees an overall relative phase 
shift in time (from the ground state energy of the zeroeth 
state), a slight dipole oscillation (from occupation of the 
first excited state) and breathing (from the second). Af- 
ter one oscillator period (310 ms in the case shown), ^2 
replicates its original value at the switch-off. 

In this case, the dynamics are quite gentle, and so the 
spatial scales of the phase and density features are always 
quite large compared with the absorption length . 
The fidelity of the writing and output of the information 
on the probe field (shown in Figs. UJd-e)) is correspond- 
ingly better than in the previous example. Additionally, 
because the |2) is trapped near the BEC center, this case 
avoids problems associated with ip2 occupying the region 
z < z c . 

Because of both the ease of analyzing the evolution 
and the high fidelity of outputting the information, we 
expect this case to be well-suited to controlled processing 
of optical information. For example, if the input pulse 
created a wavefunction ip2 corresponding to the ground 
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FIG. 6: Strong probe case: Formation of vector soli- 
tons. Dynamics resulting from using a relatively strong 
probe input fl p0 = (2tt) 5.7 MHz, Q c0 = (2tt) 8 MHz. 
In this case we have JV C = 0.5 x 10 6 Rb-87 atoms in a 
uj z = (2tt) 21 Hz trap and use /13 = f 2 z = 1/12, and 
an = a,22 = 5.36 nm corresponding to the system |1) = 
\5S 1/2 ,F = 1,M F = -1), |2> = \5S 1/2 ,F = 1, Mf = +1), 
and 1 3) = |5P 1/2 ,F = 2, Mf = 0). We artificially set 
a\2 = 1.04 an to exaggerate the phase separation dynam- 
ics, (a) Evolution of density in |2) N c \ip 2 (z, t)\ 2 /A shows the 
development and interaction of two vector solitons. The dot- 
ted lines indicate the times plotted in (c). (b) The total 
density N C (\M^, *)| 2 + \ik{z,t)\*)/A = N^l/A remains al- 
most constant in time, (c) The densities N c \ipi(z,t)\ 2 /A and 
N c \Mz,t)\ 2 /A (now on the same scale) are plotted as thin 
and thick solid curves, respectively, at the times indicated. 
The phases 0i, 02 are plotted as dotted and dashed curves. 



state of the oscillator potential, the evolution would re- 
sult only in a homogenous phase shift, proportional to 
the storage time, allowing long storage of the information 
or introduction of controllable phase shifts. By choosing 
the pulse lengths differently so several oscillator states 
are occupied, one could achieve linear processing or pulse 
reshaping. 

Furthermore, note that in this example there is a small 
but discernable dip in the |1) density at 51 ms, indicating 
some nonlinearity in the evolution of ip2 ■ One can tune 
this nonlinearity by varying the probe to coupling ratio, 
leading to nonlinear processing of the information. 
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D. Strong probe case: Two-component solitons 

The two-component dynamics are even richer when one 
uses strong probe pulses so nonlinear effects become very 
evident in the evolution. In such a case, the qualita- 
tive features of the dynamics will be strongly effected 
by whether or not the relative scattering lengths are in 
a phase separating regime an > V a n a 22 |3ll 132^. Ex- 
periments have confirmed that the |1), |2) studied in our 
previous examples in Rb-87 are very slightly in the phase 
separating regime |16| . Relative scattering lengths can 
also be tuned via Feshbach resonances ■ 

Figure [B] shows the evolution of the two components 
following the input and stopping of a stronger probe 
(f2 p o = 0.71 57 c o). In this case we chose our levels to 
be |1) = |55 1/2 ,F = 1,M F = -1), |2) = \5S 1/2 ,F = 
1,M F = +1), and 1 3) = |5P 1/2 ,F = 2,M F = 0), which 
would require an optical trap J2J] in order to trap the 
|1) and 1 2) equally. This level scheme is advantageous in 
the strong probe case since /13 = / 2 3 which, as we will 
discuss in more detail in Section IV CI improves the fi- 
delity of the writing and output processes. We chose the 
scattering length an = 1.04 an, higher than the actual 
background scattering length, to exaggerate the phase 
separation dynamics. One sees in Fig. |Bfa),(c) that over 
a 30 ms timescale, the phase separation causes the den- 
sity in 1 2) to become highly localized and dense. This 
occurs because the scattering length ai 2 > an causes |1) 
atoms to be repelled from the region occupied by the |2) 
atoms, and in turn the |2) atoms find it favorable to oc- 
cupy the resulting "well" in the |1) density. These two 
processes enhance each other until they are balanced by 
the cost of the kinetic energy associated with the increas- 
ingly large spatial derivatives and we see the formation 
of two-component (vector) solitons 01 • In the case here, 
two solitons form and propagate around the BEC, even 
interacting with each other. The alternating grey and 
white regions along each strip in Fig. Efa) indicate that 
the solitons are undergoing breathing motion on top of 
motion of their centers of mass. Fig.HJb) shows that the 
total density profile ipQ varies very little in time. It is the 
relative densities of the two components that accounts 
for nearly all the dynamics. 

We found solitons formed even in only slightly phase 
separating regimes (an/ ' \Ja\\CLT2 > 1.02). The number 
of solitons formed, the speed of their formation, and their 
width were highly dependent on this ratio as well as the 
number of atoms in |2). The ability of stopped light 
pulses to create very localized two-component structures 
seems to be a very effective method for inducing the for- 
mation of vector solitons, which has hitherto been un- 
observed in atomic BECs. A full exploration of these 
dynamics is beyond our scope here. 

These non-trivial features can also be written onto 
probe pulses as shown in Fig. [3a). Note that the condi- 
tion for coherent revivals laA^ <§; a a does not depend 
on the weak probe limit. Therefore, the writing process 
is still primarily a coherent process, and (|12|) is still well 




FIG. 7: Writing solitons onto probe field, (a) Writing 
of vector solitons onto probe and coupling fields directly after 
a switch-on after a storage time of r st = 110 ms and (b) probe 
output. The plots use the same conventions as Fig. [^d-e). 

satisfied. However, we see in Fig. Ufa), the coupling field 
is strongly affected by the writing and is far from ho- 
mogenous, meaning the weak probe result l|13|) can not 
be used. Nonetheless, we see that the qualitative features 
of both the density and phase of ip2 have still been trans- 
ferred onto Jlp. The writing in the strong probe regime 
will be studied in more detail in Section Ivl 

In Fig. Efb) we plot the resulting output probe pulse 
in this case. As in the weak probe case there is some 
degradation in the subsequent propagation due to atten- 
uation of high frequency components. In addition, unlike 
the weak probe case, there are nonlinearities in the pulse 
propagation itself which causes some additional distor- 
tion. Section IY CI also addresses this issue. Even so, we 
again see a very clear signature of the two solitons in the 
both the intensity and phase of the output probe. 

We have performed some preliminary simulations in 
traps with weaker transverse confinement and which cal- 
culate the evolution in the transverse degrees of free- 
dom. They show that the solitons can break up into 
two-component vortex patterns via the snake instability 

SHI- 

We chose our parameters to be in the phase sepa- 
rating regime here. Just as in the weak probe case, 
other regimes will lead instead to much gentler dynam- 
ics. For example, it has been shown |35| that when 
«12 < an = Q-22 then there exist "breathe-together" 
solutions of the BEC whereby complete overlap of the 
wavefunctions V'IjV^ persists. In this case the nonlin- 
ear atomic interaction can lead to spin-squeezing |37l ] . It 
would be extremely interesting to investigate a probe- 
revival experiment in such a case, to see to what extent 
the squeezed statistics are written into probe pulse, pro- 
ducing squeezed light |38j |. This would require taking 
both the light propagation and atomic dynamics in our 
formalism beyond the mean field. Steps in this direc- 
tion have been taken in [y, |39j, but these analyses are 
restricted to the weak probe case. 



V. QUANTITATIVE STUDY OF WRITING 
AND OUTPUT 

The above examples demonstrate that a rich variety 
of two component BEC dynamics can occur, depending 
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on the relative interaction strengths and densities of the 
states involved. They also show that remarkably compli- 
cated spatial features in both the density and phase can 
be written into temporal features of a probe light pulse 
and output. The purpose of this section is to quantify 
the fidelity with which wavefunctions can be written onto 
the probe field. 

To do this we consider an example of a two-component 
BEC with a Gaussian shaped feature in ip%, with pa- 
rameters characterizing length scales and amplitudes of 
density and phase variations. Note these Gaussian pulse 
shapes are relevant to cases in which one might perform 
controlled processing, e.g. Fig. Switching on a cou- 
pling field then generates and outputs a probe pulse with 
these density and phase features. We will calculate this 
output, varying the parameters over a wide range. In 
each case, we will compare the output to what one ex- 
pect from an "ideal" output (without dissipation or dis- 
tortion), and calculate an error which characterizes how 
much they differ. Using the analysis of the switching pro- 
cess in Section Hill and the USL propagation in Section ITT1 
we also obtain analytic estimates of this error. 

In the weak probe case, we find a simple relationship 
between ip2(z, t on ) and the output fip OUt ' (t). We find that 
to optimize the fidelity in this case one should choose 
the Gaussian pulse length between two important length 
scales. The shorter length scale is determined by EIT 
bandwidth considerations, which primarily contribute er- 
ror during the output. On the other hand, for large pulses 
comparable to the total condensate size, the error during 
the writing process dominates due to ^2 partially occu- 
pying the condensate entering edge region (z < z c ). 

In the strong probe case, we find that the relation- 
ship between ip2(z, t on ) and Op OUt ^ (t) is more complicated, 
however, when /13 = /23, one can still make an accurate 
correspondence. We find the fidelity in this case depends 
only weakly on the probe strength. Conversely, when 
/13 7^ /23, even a small nonlinearity causes additional 
phase shifts and distortions, making this correspondence 
much more difficult, leading to higher errors for stronger 
probes. For this reason, a system that uses, for example 
|1> = \F = l,M F = -1), |2> = \F = l,M F = +1) (as 
in Fig. UJ would be preferable when one is interested in 
applications where there are comparable densities in the 
two states. 

By considering the switch-on and output processes, 
starting with an arbitrary initial two-component BEC, 
we emphasize that this method of outputting the atomic 
field information onto light fields works regardless of how 
the BEC state was generated. One could prepare a BEC 
in a coherent superposition of |1) and \2) by any avail- 
able method including, but not limited to, inputting a 
slow light pulse. 



A. Deriving tp2 from the written and output probe 

field. 

Assume we have a two component B EC with wave - 
functions ipi , "02 (and again define ipQ = y/\ipi\ 2 + IV^I 2 ) ■ 
For now we assume /13 = fa (so an = 0, see Eq. @), 
as we treat the /13 7^ /23 case in Section TV CI Upon a 
rapid switch-on of a coupling field with amplitude fi c o, 
inverting 10 shows the pattern written onto the probe 
field will be: 

fiM=-^ D + ^ (16) 

In the ideal limit, where the spatial variations of the 
wavefunction are sufficiently small, o.na —> and the 
second term vanishes. Furthermore, lllfl shows in this 
case that the dark field intensity is constant \£Id\ = Q c o- 
So in the ideal limit we have for the amplitude: 

| n (rev-ideal)| = HMft c0 (1 7 ) 

y Wo 

In practice this will be a good approximation so long as 
the inequality |ckjva 1 1 ^1 1 ^ c^aIV^I is satisfied. The ex- 
tent to which the second term in cannot be neglected 
will determine the error between the ideal Qp CV ldcal ) and 

actual flp ev ^ output. Note that in all the cases studied 
in Section llVl ipo{ z , t) was nearly constant in time, being 
always well approximated by the original ground state 
tp\ (z). Thus ipo can be considered a known function of 
z. 

Turning now to the phase, the simple relationship 
02 — 01 — 0p ~ 4>c + 7T is always satisfied in the dark 
state (see (|12|) 'l. If C was just a constant then the probe 
phase would simply reflect the relative phase of the two 
wavefunctions P = 02 — tpi + n (choosing our phase 
conventions so <j) c = 0). This is indeed the case in the 
weak probe limit, as then the coupling field phase <f> c is 
unaffected by the atomic fields upon the switch-on. How- 
ever, this is not necessarily so in the strong probe regime. 
To calculate the true phase shift, we calculate phase of 
flo from ijTTjl. In the ideal limit <I>d{z) — f dz'cti(z') 
(choosing our phase conventions such that cf>i{zi n ) = 0). 
Defining 021 = 02 — ^i and using the definition of a; © 
and l|17[) . we find the relationship: 

4, 21 (z) = 4 rov ~ idcal '(z) + ^-4 nl) ; 

(18) 

The equation is written in this way (the atomic phase 
in terms of the probe phase and intensity) because our 
purpose is to derive the phase pattern 02 1 based on the 
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observed probe field, rather than vice versa. The non- 



linear correction 



is fact the phase imprinted on the 



coupling field <f> c during the switch-on. 

However, in an experiment, we do not have direct ac- 
cess to the revived probe dp v \z). Rather, we observe 
the temporal output Ct p out \t). In practice, the intensity 
of this quantity can be measured with a PMT, while the 
phase pattern could be measured by beating it with a 
reference probe field which did not propagate though the 
BEC. In the ideal limit (the absence of attenuation or 
distortion) the relationship (|14H will link the observed 



output to the revived probe field £1 
the "ideal" output via 



(rev] 



. Thus we define 



^(rev-idcal)^ 



Q(out- 
P 



idcal) (^out)-r(z)+t on ) (19) 



In practice, additional absorption events and distortion 
can occur during the output. If variations of the relative 
amplitude ^2/^1 have some characteristic scale then 
time features with a scale r ~ L^/V g will be introduced 
into the probe Cl p . Using our result for the bandwidth 
induced attenuation Q yields an estimate for transmis- 
sion energy (that is, the energy of the actual output pulse 



flp OUt> relative to an unattenuated output of Cl p ' 

T 

(3 



(rev) \ 



DjZont) ~ D(Zp) 

(L^a A (z p )) 2 



(20) 



For this simple estimate we evaluate oia and D at the 
location of the center of the pulse z p , though one could 
also construct more sophisticated estimates by integrat- 
ing over the spatial distribution of the pulse. Just as in 
our discussion below Eq. JSJ, the temporal width is in- 
creased by a factor T _1 and the peak intensity reduced 
by T 2 during the propagation out. 

We now apply these findings to a Gaussian pulse, with 
ip2 assumed to be of the form, 



ipziz) = ipo(z)A 2 exp[ - —J exp 



1L\ 



.A, 



-erf 



(21) 

Note tp2 has a density feature of amplitude A\ and length 
L2 and a phase feature of amplitude A^ 2 and length L^ 2 . 
For ipi (z) we choose the amplitude so the total density 

matches the ground state for the trap 4>o = IV'i G ' > | 2 and 
in the phase we put in a shift with some amplitude Afa 
and length L^ >1 . An example is shown in Fig.|S{a). These 
parameters will be varied throughout this section to learn 
how they effect the writing and output. We choose &ui z = 
(2tt) 20 Hz trap (A = tt(5 ^m) 2 ) with N c = 2.0 x 10 6 Rb- 
87 atoms. We use /13 = /23 = 1/12, corresponding to the 
system |1> = \5S 1/2 ,F= 1,M F = -1), |2> = \5S 1/2 ,F = 
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FIG. 8: Writing and outputting atomic field informa- 
tion onto the probe fields, (a) A test case with a den- 
sity feature of length L2 = 10 /im of amplitude A2 = 0.5 
and phase features of A ( j >2 = —0.2n, Lfa = 5 /xm (see 
Eq. 1211 1 and A$ x — 0.5 7r, L$ 1 — 5 ^m. The conven- 
tions for the curves are the same as Fig. OJc). (b) Upon 
switch-on of the coupling field with f2 c o = (2tv) 8 MHz, 
we show the dark intensity \Q.d(z)/0. c o\ 2 (thick solid curve) 
and coupling intensity |fl c (z)/f2 c o| 2 (thin solid). The dot- 
ted, dashed, and dot-dashed curves show, respectively, the 
phase profiles <j> c , 4>d and $21 . (c) The normalized probe in- 
tensity profile |f2p V | 2 /fi 2 (solid curve) is indistinguishable 
from the expected profile |f2 p rcv_ldca1 ' | 2 /f2 2 o (thick dashed) 
(see Eq. 1171 1. The thin solid curve shows 4> P ICV ^ and the 



dot-dashed curve shows fei + it, while the thin dashed curve 
shows the expected probe phase profile with the nonlinear 
correction (j>p ev ldca1 ' = 21 4- <^p n1 ^ + ty 11811 . For reference, 
the dotted curve shows the total density i/jq profile, (d) The 
temporal pattern of the output probe intensity If^ "*- 1 1 2 /£} 2 
(thick solid curve) and phase 4> p out ^ (thin solid), versus an 
estimate based on perfect propagation of the revived pulse 
f2 p rcv ' according to 11411 (dashed and dotted curves). The 
dot-dashed curve shows the ideal output attenuated by T 2 
(see (I20H and the discussion afterward) using the estimate 



-1/2 



1,M F = +1), and |3> = \5P 1/2 ,F = 2,M F = 0) in Rb- 
87. 

In the example of Fig. 02a), A 2 = 0.5 is not partic- 
ularly small, so the weak probe limit can not be as- 
sumed. Figure 02b) shows the spatial profiles of the 
dark field fljj and coupling field fl c immediately fol- 
lowing a fast switch-on (r s = 0.1 fis). We see very 
little attenuation of flu occurs across the BEC, as 
a N A ~ 0(A 2 L2,iA 2 A t f )2 L t j }2 ,iA2A t f )1 L r f, 1 ) is quite small 
compared with a a- Translating this back into the Q c , fl p 
basis, predicts that f2 c acquires a dip in intensity with 
a height proportional to the density \ip 2 \ 2 ■ This behavior 
is indeed seen in the figure. 

The phase difference 4>2i is plotted as the dot-dashed 
curve in Figure |SJb). One sees that in the region where 
<f>2i is inhomogeneous, a small phase shift, equal to (j)^ 
l|18p. is introduced in the coupling field (dotted curve). 
Again this shift only arises in the strong probe regime. 
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The phase shift in the dark field <po is plotted as the 
dashed curve. 

Figure H[c) then compares written probe field to our 
ideal limit predictions 1|17JI - 118|) in this example. The nu- 
merically calculated intensity |f2p rcv ' ) | 2 (solid curve) is al- 
most indistinguishable from our prediction |f2 p ICV ldeal )|2 
(|17fl (dashed curve). The phase written onto the probe 

field </>p rov ^ is also very close to the ideal limit prediction 
(|18fl . One sees that including the nonlinear correction 

</>p n1 ^ to the simpler estimate <p 2 i + tt (dot-dashed curve) 
is important in making the comparison good. 

Figure IHtd) then shows intensity (thick solid curve) 
and phase (thin solid) of the output pulse Q p out \t). For 
comparison the dashed and dotted curves show, respec- 
tively, the output expected from an unattenuated and 
undistorted transmission of the revived pulse fl p rcv \z) 
via l(T§l) . One sees a very good agreement in the phase 
pattern, while there is a visible reduction in the intensity, 
due to bandwidth considerations. The dot-dashed curve 
shows the our estimate with the estimated reduction in 
intensity T 2 (see (J^OJ) using a characteristic length scale 
L^f, calculated from a quadrature sum of contributions 
from the amplitude and phase features in (|21ll (the ex- 
pression is given in the caption). 

We have thus calculated a method by which the output 
probe pulse can be solely used to calculate the relative 
density and phase of the wavefunctions which generated 
it, and demonstrated the method with a generic example. 
Furthermore, we have identified the leading order terms 
which will cause errors in these predictions. In particular 
we have seen in our example in Fig. IHtd) that including 
the expected bandwidth attenuation accounts for most 
of the deviation between our ideal predictions and the 
actual output pulse. 



B. Quantifying and estimating the fidelity. 

We now quantify the deviations from our predictions 
(|17I19|) for our example i|21|) , varying the length and am- 
plitude parameters over a wide range. These results can 
be directly applied to pulses which are approximately 
Gaussian (as in Fig. EJ- The results here should also pro- 
vide a good guide to the expected fidelity in more com- 
plicated cases so long as the length scale and amplitude 
of features can be reasonably estimated. 

To quantify the deviation from our ideal case predic- 
tion i|17|) - (|18fl we define the write error: 







2 


JzT dZ ' 


^(rev) _ ^ (rev-ideal) | e i(0 21 +0Cnl) +7r ) 





/ Zout dz'\n { p Icv) \ 2 

(22) 

We plot this quantity in a series of cases with different L 2 
in Fig. EI a) (circles) in a case with amplitude A 2 = 0.5 



and no phase profiles. This is compared with a calculated 
prediction for E w (solid curve) based on p6|l where we 
calculate VLa with 111 Oil and calculate the small attenua- 
tion of with (|llfl . These errors grow as L 2 becomes 
comparable to = 0.2 /im, according to the discus- 
sion in Section IIII Al We see the agreement between 
the analytic and numerical estimates is quite good for 
small to moderate L 2 , confirming that this is the leading 
source of error. However, when L 2 becomes compara- 
ble to the total BEC size (the Thomas-Fermi radius is 
R z = 44 /im here), we begin to sec additional errors be- 
cause ip2 becomes non-zero at the BEC edge. Thus, we 
see that for a given BEC length and density, one must 
choose a^ 1 <C L 2 *C R z to minimize the write error E w . 
Because this edge effect depends on the density structure 
near the BEC edge, it is difficult to estimate analytically, 
but it is related to the fraction of the pulse in or near the 
region z < z c . For pulses near the BEC center, this will 
depend on L 2 /R z , whereas for pulses far off center, it will 
also depend on the location of the center of the pulse. 

We also plot a series of cases keeping L 2 = 10 /xm and 
A 2 = 0.1 but varying the length scale of a phase 
profile of amplitude — 0.75 n. One sees similar but 
slightly smaller errors E w in this case. 

We define the error accumulated as flp^ propagates 
out using l|14|) : 

g ^gV) - 4° Ut) (r(z out ) - r(z'))\ 2 

dz'i4 rcv) (zoi 2 

(23) 

In Fig. E{b) we plot this quantity for the cases corre- 
sponding to Fig.|5fa). For comparison, for the series with 
no phase shift, we calculated an estimate based on the ex- 
pected attenuation T 2 and spreading T _1 of a Gaussian 
pulse, due bandwidth considerations (see Eq. ifflTjl and 
subsequent discussion), choosing L^p — L 2 . In the limit 
of small j3 this calculation yields E ont ss (5/16)/3 2 . This 
is plotted as a solid curve and we see good agreement 
with the numerical data. 

The error E out is seen to dominate E w for small L 2 , 
due to the fact that the large optical density effects the 
former (D(z out ) = 617 in the case plotted). Thus our an- 
alytic estimate E ou t — (5/16)/? 2 is a good estimate of the 
total error. But for larger L 2 the edge effect in E w be- 
comes important. To minimize the total error E w + E out 
we should choose L 2 so the error E out — (5/16)/3 2 is 
comparable to the edge effect error. In Fig. the op- 
timal length is L 2 ° pt ' ) w 17.5 [im and the total error 

is .Bout + E w « 0.0009. The scaling of L 2 ° pt) with the 
condensate size R z is difficult to estimate. Assuming 
it roughly increases as L 2 ° pt ' ) oc R z , then l|20l) shows 
/3(°pt) oc l/D(z out ), giving us a guide as to the improve- 
ment in fidelity we can expect by increasing the total 
optical density of the condensate. 

The insets show the plots of the ideal output of 
flp CV \z) (solid curves) versus actual outputs Q p out \t) 
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FIG. 9: Writing and outputting atomic field informa- 
tion onto light fields, (a) Circles show the write error 
e w dm versus the length L2 for BECs with no phase jumps 
(j4 02 = Afa = 0) and with A 2 = 0.5 (see Eq. (J2TJ)- The 
solid curve shows an analytic result based on the analytic 
expressions 116L (1101 and 1111 . The squares show a case 
holding L2 = 10 /im, A2 = 0.1 and introducing a phase fea- 
ture A,f, 2 = 0.75 7r with varying length scale L</, 2 . (b) The 
output error -B ou t 1231 for the same cases as (a). The solid 
curves showing the prediction (5/16)/3 2 discussed in the text 
(choosing = L2) for the case with no phase jumps (cir- 
cles). The insets show the output in the cases indicated with 
arrows, with the same plotting conventions as Fig.|HJd) (in the 
lower inset on the right <j> p is homogenous and not plotted), 
we again use u x = (2n) 20 Hz, A = 7r(5 ^m) 2 , iV c = 2.0 x 10 6 
and /13 = /23 = 1/12, corresponding to the system |1) = 



|55 1/2 ,F = 1,M F = -1), |2) = \5S 1/2 ,F 
and |3) = \5P 1/2 ,F = 2, M F = 0> in Rb-87. 



l,Mj 



+1), 



(dotted curves) in two of the cases. In one case, with no 
phase jumps and a small L 2 , we see an overall reduction 
in amplitude and slight spreading. In the other, with a 
larger L 2 but a phase shift with a small length scale , 
we see that the attenuation is localized in the middle of 
the pulse. This is because it is the components which con- 
tain the sharpest phase profile, near the middle, which 
are most severely attenuated. Generally, in cases with 
complicated spatial features, one must be aware of this 
potential for local attenuation of the sharpest features. 
Note, however, that in many cases of interest, such as the 
interference fringes (Fig. 0} and solitons (Fig. HJJ), large 
phase shifts occur primarily in regions of low density, 
meaning these features often survive during the output 
propagation. 

Our analysis of switching process in Section IIIII ap- 
plies even with a large amplitude A 2 , as was seen in the 
A 2 = 0.5 case in Fig. |S] In Fig. ITDT a) we plot the results 
of a series of simulations with A 2 varying all the way up 
to A = 0.9 and see only a very small impact on E w for 
A2 < 0.8. Most of our results for the propagation of USL 
pulses, e.g. Q and (J5J, rely on the weak probe limit. 
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FIG. 10: Writing and outputting atomic field informa- 
tion onto light fields, (a) The write error E w \22l versus 
A-2 for the case L 2 = 10 fim, L ( p 1 = L^, 2 = 5 fim, A <j>l = 
0.57T, Afo = 0.27T. (b) The output error i5 ou t ll2lfl l in the 

same cases. The solid line shows the estimate (5/16)/3 2 (see 
Eq. OOJ) using Lnp = (L2 2 + A^L^ + A^Lg)- 1 ' 2 . The 
condensate parameters are the same as in Fig. [§] 



Surprisingly though Fig JlOf b) shows that the output er- 
ror E out is independent of A 2 through A 2 = 0.9. 



C. Effect of unequal oscillator strengths 

The lack of distortion in the strong probe cases dis- 
cussed so far in this section is a result of us choosing a sys- 
tem with equal oscillator strengths /13 — f 2 3- The impor- 
tance of the relative size of the oscillator strengths can be 
understood by considering the adiabatons (the changes 
in the output coupling intensity seen in Fig.[5fa,e)). The 
adiabatons arise because of the coherent flow of photons 
between the two light fields and their amplitude is deter- 
mined by the requirement that the sum of the number of 
photons in both fields is constant upon such an exchange. 
In a case where fL ~ fL one must take into account that 



it is the total field strength 



\n p \ 2 



|S! r | 2 , rather 



than simply the input coupling strength Qj, Q which en- 
ters the numerator in our equation for the group velocity 
V g (z) (0J. However, the expressions for the Rabi fre- 
quencies involve the oscillator strengths of the relevant 
transitions. Thus, preserving the total number of pho- 
tons at each point in space will result in a homogenous 
|£!| 2 = f^ if and only if /13 = f 2 $. During the time that 
the probe is completely contained in the BEC, |f2| 2 will 
differ from f2 2 in the vicinity of the pulse according to: 



\n(z,t)\ 2 = n 2 c0 + \n p (z,t)\ 2 1-^ 



Jj3 
fl3 



(24) 



This leads to a nonlinear distortion, whereby the cen- 
ter of the pulse (where the probe intensity is greatest) 
tends to move faster or slower than the front and back 
edges, depending on the sign of term in parentheses. 
This behavior is quite clearly seen in Fig. 1111 where 
we compare a delayed weak probe regular USL output 
pulse to cases with strong probes 51 p o = ^ c o and vary 
f 23 = 1/3,1/2,2/3 (/13 = 1/2 in all cases). Each time 
we chose the coupling intensity so that that the Rabi 
frequency is J7 cn = (2tt) 8 MHz. One sees that when 
/13 = /23 there is no asymmetrical distortion, though we 



17 



1 




FIG. 11: USL with a strong probe. The dotted curve 
shows a strong probe pulse amplitude with peak value fl p o = 
fl c o = (2ir) 8 MHz input into a sodium cloud with an optical 
density D(z ou t) = 407. We use an oscillator strength /13 = 
1/2 but vary the coupling oscillator strength f23. The output 
f2p OUt ' /fi p o is shown in the cases /a3 = 1/2 (thick solid curve), 
/23 = 2/3 (dot-dashed curve), and /23 = 1/3 (dashed curve). 
For reference, the thin solid curve shows output for a weak 
input probe il p0 = (2n) 1.4 MHz. 



get a slight spreading and a shorter delay time. However, 
in the other two cases, there is asymmetry in the output 
pulses, consistent with the picture that the center of the 
pulse will travel slower or faster than the edges according 
to (|2"4"|) and I0J (with Q^o replaced by |f2| 2 in the numera- 
tor). The result is a breakdown of our estimate (|19fl . As 
this is a propagation effect, the absolute magnitude of 
this effect should roughly be linear with optical density 
and the relative distortion will go as the optical density 
divided by the pulse length. This is the primary reason 
it would be advantageous to use the levels we chose in 
our simulation in Fig. [71 if there is a strong probe. 

We just outlined the reason that the equal oscillator 
strengths can be important in reducing distortion during 
propagation and thus E out . It turns out that the writing 
process is also more robust when /13 = / 2 3. The rea- 
son is the presence of the additional term au JHJl. This 
term is quite small in the weak probe limit, but in the 
strong probe case leads to an additional phase shift, not 
accounted for in i|18|) . which depends in detail on the 
pulse amplitude and structure. 

Figs. I12f a-d) shows a case, similar to Fig. [5J but with 
/13 = 1/12, / 2 3 = 1/4 (as in the cases in Figs. I4l5|l . and 
with a fairly small nonlinearity A2 = 0.25. One sees in 
Fig. I12f b) that there is a dip in Hp in the region of the 
probe due to the a±2 term. Fig. I12f c) shows the written 
probe pulse. There is a small but discernable difference 
between the predicted </>p rev_ldcal - ) (|18l) and actual <pp ev ^ 
phase. Fig. I12f d) then shows the asymmetric distortion 
which develops as the strong probe propagates out. One 
sees the phase jump is distorted in addition to the ampli- 
tude. Figs. ll2f e-f) demonstrate how these effects lead to 
substantially higher errors E w and E oni as A2 becomes 
larger. From our estimate from the previous section we 
expect an output error of E out — 0.0147 and we see in 
Fig. I12f f) that the distortion effect leads to errors higher 
than this when A2 > 0.3. 

We thus conclude that fidelity of both the writing and 



a b 





A2 A2 

FIG. 12: Writing information with unequal oscillator 
strengths, (a-d) Plots are for the same parameters as Fig.|H| 
except with /12 = 1/12, / 23 = 1/4 and A 2 = 0.25. (e-f) 
Write E w and output E out errors versus A2 for the parameters 
1/2 = 10 fim, L ( / >1 = L^ 2 — 5 /jm, A ( p 1 = 0.5 tv, A^, 2 — — 0.2 7r. 

output is compromised in a system with unequal oscil- 
lator strengths. However, this only comes into effect in 
cases where the nonlinearity is important. In weak probe 
cases neither of these effects is important and the fidelity 
can still be estimated with the analysis of the /13 = /2s 
case in Sections IV AIV Bl 

VI. OUTLOOK 

In conclusion, we have established first several new re- 
sults regarding the fast switching of the coupling field 
in light storage experiments. We found the switching 
could be done arbitrarily fast without inducing absorp- 
tions so long as the probe is completely contained in the 
atomic medium (and V g <C c which is usually the case 
in practice). We have also seen that when the switching 
is slow compared to the excited state lifetime (~ 25 ns), 
the probe smoothly follows the temporal switching of the 
coupling field, while in the other limit the probe am- 
plitude undergoes oscillations which damp out with this 
time scale (see Fig.[3Jl. 

Next, we saw that these stopped pulses can induce 
novel and rich two-component BECs dynamics during 
the storage time. Both the relative scattering lengths 
of the states used and the probe to coupling intensity 
ratio have a strong effect on the qualitative features of 
the dynamics. In the weak probe case, the |2) con- 
densate sees en either an effective repulsive hill (when 
ai2 > an, see FigQJb)) or harmonic oscillator poten- 
tial (when ai2 < an, see FiglSJa)). The characteristic 
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timescale for the dynamics is given my the chemical po- 
tential of the BEC (in our cases ~ (2ir) 1 kHz) times the 
relative scattering length difference \a\ 2 — a\\\/ax\. In 
the latter case, the resulting evolution can be easily calcu- 
lated by decomposing the pulse into the various harmonic 
oscillator eigenstates. Thus it is possible to choose input 
pulses which preserve their density over time or undergo 
predictable reshaping such as dipole sloshing or breath- 
ing, allowing controlled storage and linear processing of 
optical information. Inputting stronger probes will add 
nonlinearity to the evolution, making nonlinear process- 
ing possible. Very strong probes and phase separating 
scattering lengths a\ 2 > ^11022 lead to the formation and 
motion of vector solitons, which have not been observed 
in BECs to date. 

We then showed that switching the coupling field on 
after the dynamics writes the various density and phase 
features of the wavefunctions onto revived probe pulses. 
This was seen qualitatively for various examples (Figs.QJ- 
|7J) . A precise relationship between the wavefunctions and 
output pulses was found (jl7|l - 119|l . meaning the output 
probe pulse can be used as a diagnostic of the relative 
density and phase in the BEC wavefunctions. We have 
also identified sources of attenuation and distortion in 
the writing and output processes and quantitative errors 
were calculated for a wide range of length and amplitudes 
of Gaussian pulses (Figs. I9I1U|) . As we saw there, the er- 
ror during the output dominates for shorter pulses and 
is (5/16)/3 2 (see Eq. (j20J))- For longer pulses effects due 
to the condensate edge introduce errors into the writing 
process. Balancing these two considerations one can op- 
timize the fidelity, which improves with optical density. 

For the strong probe case, we found that for equal os- 
cillator strengths (/13 = 723), one could still relate the 
wavefunctions to the output pulses by taking into account 
an additional nonlinear phase shift l|18fl . The fidelity of 
transfer of information was virtually independent of the 
probe strength Sl 2 even when Q p0 ~ Q c o (see Fig. 1117)1 . 
For unequal oscillator strengths (/13 7^ 723) we found 
strong probes lead to additional features in the phase 
pattern during the writing process and in distortion dur- 
ing the propagation during the output, leading to much 
higher errors (see Fig. I12fl . 

Looking towards future work, we note that it is triv- 
ial to extend the analysis here to cases where there are 
two or more spatially distinct BECs present. So long as 
they are optically connected, information contained in 
the form of excitations of one BEC could be output with 
this method and then input to another nearby BEC, lead- 
ing to a network. In the near future, there is the exciting 
possibility of extending these results beyond the mean 
field to learn how using this method of writing onto light 
pulses could be used as a diagnostic of quantum evolu- 
tion in BECs. A particular example of interest would 
be to investigate the spin squeezing due to atom-atom 
interactions in a two-component system and the subse- 
quent writing of the squeezed statistics onto the output 
probe pulses. Furthermore, we expect that performing 



revival experiments after long times could be used as a 
sensitive probe of decoherence in BEC dynamics, similar 
to the proposal in [Sof. Lastly, there is the prospect of 
inputting two or more pulses in a BEC and using control- 
lable nonlinear processing from atom-atom interactions 
to design of multiple bit gates, such as conditional phase 
gates. We anticipate the results presented here can be 
applied to these problems as well as to applications in 
quantum information storage, which require the ability 
to transfer coherent information between light and atom 
fields. 



APPENDIX: ADIABATIC ELIMINATION OF ip 3 
AND FASTER SWITCHING 



We arrived at Eqs. Q-@ by first considering a system 
of equations with all three levels considered and then 
adiabatically eliminating the wavefunction for |3). For 
completeness we write the original equations here. 

The three coupled GP equations are: 



dt 



1 



^12(^2 



2 

2m dz 2 



5^3, 



r 



in 

(A.1) 

(A.2) 
(A.3) 



Note we have ignored the external dynamics of ips alto- 
gether as they will be negligible compared with T, J7 p , f2 c 
(on the order of three orders of magnitude slower). 
Maxwell's equations read: 



d , 1 d\ ./i 3 a r 



9 1 d , „ 
1 ) Q c 

dz c dt 



In adiabatically eliminating ip 3 , we assume all quantities 
vary slowly compared to excited state lifetime (~ 
16 ns in sodium). This allows us to set dips/dt — ► in 
3J and arrive at: 



^3 ~ - r= (f2 p Ipl + U c 1p 2 ) 



(A.5) 



Plugging this into !|A.l|) - (|A.2jl and !|A.4|) yields CJ-0. 
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